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ABSTRACT 

We present the discovery of an extremely bright and extended lensed source from the second Red 
Sequence Cluster Survey (RCS2). RCSGA 032727-132609 is spectroscopically confirmed as a giant 
arc and counter-image of a background galaxy at z = 1.701, strongly-lensed by the foreground galaxy 
cluster RCS2 032727-132623 at z = 0.564. The giant arc extends over ^ 38 "and has an integrated 
r-band magnitude of 19.1, making it ^ 20 times larger and ^ 3.5 times brighter than the prototypical 
lensed galaxy MS1512-cB58. This is the brightest distant lensed galaxy in the Universe known to date. 
We have collected photometry in 9 bands, ranging from u to Kg, which densely sample the rest-frame 
UV and optical light, including the age-sensitive 4000 A break. A lens model is constructed for the 
system and results in a robust total magnification of 2.04± 0.16 for the counter-image; we estimate an 
average magnification of 17.2± 1.4 for the giant arc based on the relative physical scales of the arc and 
counter-image on the sky. Fits of single-component spectral energy distribution (SED) models to the 
photometry result in a moderately young age, t — 80±40 Myr, small amounts of dust, E{B—V) < 0.11, 
and an exponentially declining star formation history with e-folding time r = 10 — 50 Myr. After 
correcting for the lensing magnification, we find a stellar mass of M* ~ 10^° M0 and a current star 
formation rate SFR< 77 Mq yr^^. Allowing for episodic star formation, an underlying old burst could 
contain up to twice the mass inferred from single-component modeling. RCSGA 032727-132609 is 
typical of the known population of star-forming galaxies near this redshift in terms of its age and 
stellar mass. Its large magnification and spatial extent provide a unique opportunity to study the 
physical properties of an individual high-redshift star-forming galaxy in great detail, opening up a 
new window to the process of galaxy evolution between z — 1.7 and our local Universe. 
Subject headings: galaxies: evolution, galaxies: high-redshift, gravitational lensing 
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1. INTRODUCTION 

Significant progress has been made in recent years 
towards the study of the formation of galaxies and 
their evolution into the population of objects we observe 
around us today. The redshift range 1.0 < z < 3.0 rep- 
resents a crucial period in this process since the peak of 
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star formation in the Universe is known to occur near 



z ^ 2 (iBlain et al.||1999l [Chapman et al.||2003l iReddy k 
Steidelli>009p — ^^ ' ' '^ — 



i'hrough the development of pre-seiection 



color criteria ([Adelberger et al. | 2004 [ [Steidel et al 
and advances 



2004 
near-lK anc 



Daddi et al. 2004a) and advances in the 
U V spectroscopic capabilities of 8-10 m class telescopes, 
growing samples of optically or near-IR selected galax- 
ies are now spectroscopically confirmed in this redshift 
range and studied extensively (see e.g. Shapley et al. 
2005; Erb et al. 2006; Kriek et al. 2008). For a complete 
understanding of the process of galaxy formation, it is 
crucial to complement these statistical results based on 
larger samples with detailed study of the stellar popula- 
tions and dynamics of individual objects. The majority 
of galaxies in the current samples are too faint for this 
purpose. The few galaxies that are bright enough to be 
studied individually represent outliers drawn from the 
extreme bright tail of the luminosity function, and are 
therefore not necessarily representative of the bulk of the 
population. The ability to study the properties of faint, 
high-redshift galaxies is one of the main science drivers 
for the construction of 30 m class telescopes. Until such 
instruments become available we can make a head start 
using gravitational lensing to increase the power of the 
current generation of telescopes. 

The lensing magnification induced by foreground 
galaxy clusters and individual galaxies has been succes- 
fuUy used to identify galaxies out to z ~ 10, opening 
up new windows into the very distant Universe ( Richard 



et al. |2008 
redshitts it 
at 1.0 < 



Bouwens et al. 20091. At more moderate 
jrings individual galaxies from the samples 



_ z < 3.0 to a flux level amenable for exten- 
sive follow-up observations at various wavelengths. The 
first and most notable example in this class is MS1 512- 
cB58. a Lyman Break Galaxy (L BG) at z = 2.73 (Yee 



mg magmnc ation ot ^ 30 ( 



et al.|[T996 Ellingson et al. 1996|), found to have a lens 



Williams & Lewis 



1996 Seitz 



et al. 1998). It has been studied extensively since its 
discovery and provides a wealth of information on the 
stellar population and dynamics of a youn g star -forming 



galaxy a. t this redshift (Pettini et al.|20d0 2002 Siana et 



al. 



2008). Other examples ot particularly bright lensed 



axies at similar redshifts include the 'Cosmic Ey e' at 



gal 
z ■■ 
'8 o'clock arc' at z 



3.07 (Small et al. 



2007 



Coppin et al. 2007), the 



T73~iAliam et al. 



2D07 



Finkel 
stein et al. 2009) and two strongly- lense d^z ^3 LBU s 
from the SDSS Giant Arcs Survey (Koester et al.|[20T0 ). 
These highly magnified sources currently represent the 
best places to study the individual properties of high- 
redshift star-forming galaxies in great detail. 

In this paper we present the discovery of a very bright 
and extended galaxy, RCSGA 032727-132609, spectro- 
scopically confirmed at z — 1.701 and highly magnified 
by a foreground cluster at z = 0.564 from the second 
Red Sequence Cluster Survey (RCS2; D. Gilbank et al. 
2010, in preparation). The RCS surveys were designed 
for the purpose of cluster finding via the identification 
of the linear color-magnitude relation present for early- 
type g alaxies in clusters, know n as the red-sequence tech- 
nique ( iqiadders fc Yee||2000| . The RCS2 survey has im- 
aged ~ 700 square degrees of sky in g^r and z with the 
MegaCam camera at the Canada-France-Hawaii Tele- 
scope (CFHT) on Mauna Kea in 4, 8 and 6min exposures 
respectively. Data acquisition for the survey finished in 
2008 and preliminary cluster catalogs over the entire area 
have been created and visually inspected for strong lens- 
ing signatures; the details of this search will be published 
elsewhere. The brightest and most obvious strong lensing 
system found in this search is RCSGA (RCS Giant Arc) 
032727-132609. The system consists of a counter-image 
and a giant arc extending over ~ 38 "at an Einstein- 
radius of ^ 17.8", estimated from the distance between 
the brightest knot in the arc and the brightest cluster 
galaxy. The arc has an apparent magnitude of r — 19.1, 
making it ^ 3.5 times brighter than cB58 and the bright- 
est distant lensed galaxy in the Universe known to date. 
Its large spatial extent provides unique opportunities to 
look inside a high-redshift galaxy and spatially resolve 
its substructure. 



The paper is organized as follows. § [2] presents the 
multitude of photometric and spectroscopic data we have 
assembled on this system: broadband observations in 9 
bands ranging from u to K^^ a medium-resolution opti- 
cal spectrum of the lensed galaxy to obtain the source 
redshift and redshift measurements of 49 cluster mem- 
bers to estimate the virial mass of the foreground clus- 
ter. § [3] describes an innovative method used to obtain 
accurate photometry of the source. A lens model is con- 
structed for the cluster in § |4j We fit spectral energy dis- 
tribution (SED) models to the broadband photometry of 
the source to explore its star formation history and stel- 
lar population parameters; the SED modeling procedure 



and results are discussed in §[5l RCSGA 032727-132609 
is compared to the known galaxy population at z ^ 2 in 
§ [6J Throughout this paper we will assume ^m = 0.3, 
J7a — 0.7 and Hq = 70kms~^ Mpc~^. All magnitudes 
are quoted in the AB system. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Imaging 

We have obtained imaging of RCSGA 032727-132609 
in 9 bands ranging from u to K^ on 5 different telescopes. 
Details of these observations are given below and sum- 
marized in Table [1] 

The u-band observation is part of a larger program to 
survey ~100 strong lensing galaxy clusters in order to 
identify lensed u-band dropouts and constrain the frac- 
tion of strongly lensed galaxies that lie at z> 3 (Bayliss et 
al. 2010 in preparation). We observed RCSGA 032727- 
132609 in the SDSS u-band at the 4.1m Southern As- 
trophysical Research (SOAR) Telescope in Chile with 
the SOAR Optical Imager (SOI) on 2008, October 29 in 
six exposures of 600 s each. The data were reduced and 
stacked in IRAF [J using procedures from the MSCRED 
package. 

Pre-imaging with the F0RS2 instrument at the Very 
Large Telescope (VLT) in Chile for the purpose of mask 
design for the spectroscopy described below, consisted 
of dithered integrations in the J- and B-bands totalling 
480 s and 738 s respectively. The imaging in both filters 
was acquired on 2006, August 17. 

The g- and r-band data come from the RCS2 dataset, 
as detailed above. The data were acquired in queue 
mode in semester 2005B using the prime focus Mega- 
Cam imager on CFHT. It is important to keep in mind 
that the g-filter used for MegaCam on CFHT is simi- 
lar to, but subtly different (both in central wavelength 
and width) from its namesake in the Sloan Digital Sky 
Survey (see Gilbank et al. in preparation for further de- 
tails). The processed images are used as delivered by the 
CFHT queue pipeline, although astrometric and photo- 
metric re-calibration is performed as part of the RCS2 
catalog processing. Unlike other data used in our analy- 
sis these images are single exposures with no dithering; 
cosmic rays and chip defects - neither of which are sig- 
nificant - have been identified by visual inspection and 
removed by interpolation. 

During a run in October 2009, the IMACS instrument 
and the PANIC near-IR camera on the Magellan I Baade 
Telescope were used for broadband imaging of RCSGA 
032727-132609 in the z-, J- and X^-bands. A total of 
1080 s of dithered z-band images were acquired in nine 
individual integrations with the IMACS f/2 spectrograph 
in imaging mode. The data were processed and stacked 
using standard IRAF tasks. Dithered observations in 
both J and Kg were acquired with the PANIC near-IR 
imager, totaling 2160 s in both bands. The PANIC data 
were reduced using a custom pipeline built in IRAF with 
standard techniques of dark subtration, flat-fielding, and 
sky subtraction with iterated object masking. 

Finally, an i/-band image was constructed from 7200 s 

^ IRAF (Image Reduction and Analysis Facility) is distributed 
by the National Optical Astronomy Observatories, which are op- 
erated by AURA, Inc., under cooperative agreement with the Na- 
tional Science Foundation. 



TABLE 1 

Summary of Imaging Data 



Filter 


Total Int. 


Seeing 


Date(s) 


Telescope, Instrument 




(seconds) 


(arcsec) 


Observed 




u 


3600 


0.85 


2008-10-29 


SOAR, SOI 


B 


738 


0.84 


2006-08-17 


VLT, FORS2 


9 


240 


0.74 


2005-11-04 


CFHT, MegaCam 


r 


480 


0.64 


2005-11-04 


CFHT, MegaCam 


I 


480 


0.66 


2006-08-17 


VLT, FORS2 


z 


1080 


0.73 


2009-10-21 


Magellan Baade, IMACS f/2 


J 


2160 


0.57 


2009-10-19 


Magellan Baade, PANIC 


H 


7200 


0.96 


2009-02-06 
2009-11-28 


Apache Point 3.5m, NICFPS 


Ks 


2160 


0.51 


2009-10-19 


Magellan Baade, PANIC 



of dithered observations taken on 2009, February 6 and 
November 28, with the 3.5 m telescope at Apache Point 
Observatory (APO) in New Mexico. Additional short 
exposure J- and i^Tg-band imaging was collected for cal- 
ibration purposes using the same method. We used a 
pipeline of standard IRAF tasks to dark-subtract and 
stack the dithered images. 

Figure [l] shows a composite color-image of RCSGA 
032727-132609 constructed from all 9 bands, combin- 
ing uBg, rlz and JHKs individual composites. Postage 
stamp cutouts of all 9 photometry bands are shown in 
Figure [2] The west end of the arc is complicated by two 
cluster members that fall on top of it (see § Is] and Fig- 
ure pi) ; the u-band data obviously shows this region to 
be part of the giant arc. 

2.2. VLT Spectroscopy 

The Very Large Telescope (VLT) in Chile was used to 
measure the redshift of the lensed background galaxy. 
A total integration time of 2 hours was obtained with 
F0RS2 with a multi-slit mask in service mode on the 
nights of December 27, 2006 and January 13 and 27, 
2007. The setup used the GRIS150I disperser and GG435 
filter, providing a wavelength coverage from 4450 to 
lOSOOA with a resolution of 3.45Apixel~^. The sys- 
tem was binned to 2 pixels both spatially and spectrally. 
Spectra were taken of several locations in the main arc, a 
likely counter-image, and a few possible cluster galaxies 
selected on the basis oi B — I color as measured from 
the mask design pre-imaging. Reductions were carried 
out with the ESO package esoresQand wavelength cali- 
bration was done by comparison to standard He-|-Ne+Ar 
lamp observations. Individual images were sky-line sub- 
tracted and then combined into a single image to elimi- 
nate cosmic rays. Sky-subtraction is complicated by the 
complexity of the object placement on the slits, as this 
precluded along-slit dithering. There are thus some sig- 
nificant fringing residuals at redder wavelengths. Never- 
theless, the lensed source redshift is obvious due to strong 
absorption and emission lines; these spectra also confirm 
the identity of the counter-image. The stacked ID spec- 
tra for the arc and counter- image are shown in Figure [31 
The source redshift is 1.7009 ± 0.0008, estabhshed from 
the Fell and Mgll lines in absorption and the [OII]A3727 
and the CIII1A1909 lines in emission. 



One prominent feature in the spectrum is the 
CIIIJA1909 nebular emission line, a coUisionally excited, 
semi-forbidden transition often present in local star- 
bursts. It i s found to be stronger in starbursts of lower 



metallicity (Heckman et al. 1998), most likely due to a 



decrease in the nebular electron temperature of higher 
metallicity gas, which causes more of the nebular cooling 
from coUisionally excited lines to occur in the infrared 
rather than the UV. At higher redshifts this line is usu- 
ally too weak to be identified in individual spectra, but 
it has been detecte d in composite spectra of star-forming 



galaxies at z ~ 2 (Erb et al. 2006 
an d in a composite LBG spectrum 



Halliday et al."2008l 
3 ([Shaplcy (M 



al 



al. [2003 ). The inverse relation between the strength ot 
ClllJ nebular emission and the metallicity of the galaxy 
remains valid at these higher redshifts. It is interesting to 
note that the GUI] line was n ot detected in the UV spec- 
trum of cB58 ( [Pettini et aIJ 2000, 2002). This emission 
line can also be used m combination with either Lya or 
CIV line intensities a s an indicator of AGN contribution 
( [Shapley et aI1[2003l ). 



^ http://www.eso.org/sci/data-processing/software/pipelines/fors observed 



2.3. Cluster Velocity Dispersion 

The cluster RCS2 032727-132623 was observed spec- 
troscopically with the GISMO instrument (Gladders et 
al. , in preparation) on the Magellan I Baade Telescope on 
January 25, 2009, using the 150 1/mm grating and the 
WB5650-9200 spectroscopic filter. A single mask cov- 
ering a roughly circular area on the sky with a radius 
'^ 1.75arcmin, which corresponds to a physical radius of 
^ 0.9 Mpc at the redshift of the cluster, was observed 
in three integrations of 2400s each (7200s total). Data 
were reduced using the COSMOS pipelincrj These data 
yielded a total of 66 reliable redshifts, ofwhich 46 are 
cluster members. Three additional cluster members have 
measured redshifts from the VLT spectroscopy discussed 
above. Figure H] shows a histogram of all measured red- 
shifts in this field. No other significant structures are 
detected in redshift space to the limits of these data; 
given the tight centering of these data on the cluster 
core this is not surprising, though it does confirm that 
no other cluster-scale halos at z < 1 contribute to the 
strong lensing along the line of sight. The cluster has 
a measured redshift of 0.5637±0.0007. The best fitting 
velocity dispersion is 988 ± 122kms~^ in the rest-fr ame, 
corresponding to a M200 ^ 1-1 x IO^^Mq cluster (Yee 
& Ellingson [2003 ), in accord with the large arc radius 
here is no evidence for substructure in the 




Fig. 1. — Composite color-image of uBg, rlz and JHKs combined imaging data, (1 X larcmin). All images were transformed to the 
pixel scale of the VLT 7-band image and the PSFs were degraded to match those of the i?-band, which has the worst seeing. North is up 
and East is left. 




Fig. 2. — Gallery of 1 x larcmin images of RCSGA 032727-132609 in u,g, B,r,I, z, J, H and Ks (top left to bottom right). All images 
were transformed to the pixel scale of the VLT /-band image and the PSFs were degraded to match those of the _ff-band, which has the 
worst seeing. North is up and East is left. 
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Fig. 3. — Normalized stacked VLT spectra of the arc (top) and 
the counter-image (bottom) at rest-frame wavelengths. Important 
spectral lines are identified. 
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Fig. 4. — The main panel shows a histogram of redshifts binned at 
a redshift interval of 0.005, illustrating the dominance of the cluster 
in redshift space. The inset shows the velocity distribution relative 
to the cluster mean redshift, binned at a velocity interval of 250 
km sec~^ (rest-frame), with the best fitting gaussian distribution 
overplotted. 



cluster in velocity space. 

^ http://obs.carnegiescience.edu/Code/cosmos 



3. PHOTOMETRY 

We have developed a self-consistent, complex IDL-code 
to obtain accurate magnitudes for extended sources in a 
series of broadband images with differing point spread 
functions (PSF). All images are first transformed to the 
orientation and scale of the VLT /-band image. An em- 
pirical, normalized PSF is created for each image based 
on a well-defined, isolated and non-saturated reference 
star. We match the PSF of each frame to the i/-band 
image (which presents the worst seeing over our sam- 
ple of 9 bands) by convolving it with an appropriate 
gaussian. This sequence of PSF-matching ensures that 
we measure the light in the same physical region of the 
source in all bands. This is a crucial requirement for 
meaningful modeling of the spectral energy distribution 
at a later stage. The PSFs of these data are sufficiently 
similar that convolution with a simple gaussian provides 
adequate matching between images. Object apertures 
are created by tracing a curve along the extended source 
and convolving it with the H-hand PSF. A series of aper- 
tures of increasing radial extent are defined as isophotes 
of this convolution. The apertures are very much non- 
circular, especially for the arc, and described by equiva- 
lent radii based on circular apertures that extend to the 
same isophotes. After a detailed sequence of sky subtrac- 
tion and outlier masking, final magnitudes are measured 
at an equivalent radius of twice the FWHM of the H- 
band image, 7.6 pixels or 1'.'92, and aperture-corrected to 
an equivalent radius of 6" based on the curve of growth 
of the PSF reference star. 



Lensed background sources are often found close to the 
center of the foreground galaxy cluster. In these gener- 
ally crowded environments, problems arise when cluster 
members lie too close to the background galaxy and fall 
inside the me asurement apertur es. We use the GAL- 
FIT package (Peng et al. 2002) to fit a Sersic profile 
to these galaxies in the reterence /-band image. This 
model is convolved and scaled as needed, and subtracted 
from each PSF-matched image. Figure [5] shows the /- 
band image before and after this model subtraction, and 
overlayed with the arc apertures used for the final mea- 
surement. In the case of RCSGA 032727-132609, we are 
presented with the unfortunate and complex situation 
where 2 cluster members fall on top of the giant arc. 
They are circled in red in the top panel of Figure [sj 
Evidently their flux contribution has to be subtracted 
from the arc magnitude, but in this case the appropriate 
scaling can not be estimated easily. We create a color 
image for each frame by subtracting a scaled version of 
the r-band to make the arc vanish. The GALFIT mod- 
els of the overlapping cluster members are then used to 
subtract any remaining positive flux at their positions 
and the original frame is restored by adding the scaled 
r-band back in. This process is unnecessary for the u- 
, B- and g-bands, were the contribution of the cluster 
members is negligible, and becomes more important for 
the longer wavelengths. We conflrm the subtraction by 
comparing the brightness of the affected part of the arc 
to the brightness of other parts (an advantage of large 
spatial coverage). Through magnitude measurements of 
visually obvious over- and undersubtracted frames, we 
estimate the uncertainty which arises from this subtrac- 
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tion: 0.01 mag for r, I and z; 0.02 mag for J; 0.03 mag 
for H; 0.04 mag for Kg. This is incorporated into the 
final photometric uncertainties. 




Fig. 5. — Top: Cutout of PSF-matched /-band; circles show 
galaxies that need to be subtracted - the red circles correspond to 
the 2 cluster members on top of the arc. Middle: Final subtracted 
/-band image. Bottom: Final /-band image with arc apertures. 



The cahbration of the images is handled carefully to en- 
sure an accurate magnitude measurement of our object 
of interest. We calibrate the g-, r- and z-bands against 
the RCS2 catalog, which itself has been calibrated to 
match the SDSS stellar locus (Gilbank et al. in prepa- 
ration). There is a 0.05 mag absolute calibration uncer- 
tainty associated with the RCS2 photometry. The B- 
and /-bands are calibra ted from the qrz-ph otom etry us- 
ing transformations by Jordi et al. ( |200 6) and 'ChonisI 
[& GaskcU ( 2008 ) . The u-band calibration follows from 
the gr-photometry via a general relation derived from 
SDSS data. The near-IR J-, H- and i^Tg-bands are cal- 
ibrated against the 2 Micron All Sky Survey (2MASS). 
The small field of view of the Magellan J and Kg im- 
ages does not include enough 2MASS stars for a ro- 
bust calibration, so we use the larger APO images to 
pick up enough 2MASS objects and link the APO zero- 
point to the Magellan zeropoint based on a number of 
fainter point sources. Each star used for zeropoint cal- 
ibration undergoes the same procedure explained above 
for its flux measurement. This eliminates any differences 
between our measurement method and the surveys we 



are calibrating against and results in the most accurate 
zeropoint determination. We add the zeropoint uncer- 
tainty and the Poisson error from the flux measurement 
in quadrature to construct final photometric uncertain- 
ties. 

Final magnitu des are corrected for galactic extinction 
(Schlegel 1998) and presented in Table [2] for the main 
arc and counter-image seperately. The photometry of the 
smaller and more regularly-shaped counter-image does 
not require the careful masking of cluster galaxies inside 
the source aperture and its measurement resembles much 
more a typical photometry measurement. After tak- 
ing into account the lensing magnifications, the counter- 
image can therefore be used as a consistency check for 
the photometry of the giant arc. The arc magnitudes 
are on average a factor of 2.34 brighter, which is consis- 
tent with the difference in their sizes due to the lensing 
magnification as explained in § [4] 

4. LENS MODEL 

The multiply-lensed system RCSGA 032727-132609 is 
composed of a highly-stretched giant arc north of the 
center of the cluster, and a much less magnified counter- 
image south of it. The separation between the giant arc 
and the counter-image is ~46". The structure of the gi- 
ant arc is slightly complicated by the lensing potential 
of two cluster-member galaxies at its west side (see Fig- 
ure p|. The giant arc is highly magnified in the tangen- 
tial airection, revealing substructure that is not evident 
in the counter-image. With the resolution of the current 
imaging data we cannot uniquely identify the emission 
knots that form this substructure as multiply-imaged 
counterparts of the same emission region in the source. 
In constructing a lens model, we are therefore forced to 
make assumptions on the identity of these features. High 
resolution data may resolve some of the discrepancy and 
allow for a more detailed model in the future. 

The lens model is comp uted using the p ublicly- 
available software LENSTDOL (jJullo et al.||2007|), in an 
iterative manner. We start by representing the cluster 
halo by a single mass in the form of a pseudo-isothermal 
ellipsoidal mass distribution (PIEMD; Limousin et al. 
2005). We use as constraints the redshift of the source, 
the coordinates of the counter- image, and a series of coor- 
dinates tracing the giant arc. The best-fit model is found 
through Markov Chain Monte Carlo (MCMC) minimiza- 
tion in the source plane. The resulting first-order model 
is not a satisfactory fit to the data, but it indicates that 
the giant arc is composed of three images of the source. 
As a second step, we break the giant arc into three differ- 
ent parts, using the observed symmetry of the substruc- 
ture and guided by the location of the critical curves of 
the initial model. In our final model, we add the lensing 
contributions of the 30 brightest cluster-member galax- 
ies, each represented by a PIEMD with parameters that 
follow the observed properties of the galaxies through 
scaling relations (see Limousin et al. 2007 for further 
description of the scaling relations). Cluster members 
are selected via their location relative to the cluster red 
sequence in a color-magnitude diagram. We let all the 
parameters of the cluster halo vary in the MCMC mini- 
mization, but introduce priors on t he v elocity dispersion 
following the observed value (see § 2.3 1. The parameters 



TABLE 2 

Final magnitudes for the arc (A) and counter-image (C). 



B 



H 



Ks 



A 
C 



19.17±0.100 
21.39±0.100 



19.14±0.070 
21.26±0.075 



19.15±0.055 
21.36±0.060 



19.07±0.065 
21.29±0.060 



19.02±0.080 
21.39±0.080 



18.90±0.065 
21.22±0.070 



18.33±0.090 
20.72±0.080 



18.36±0.075 
20.83±0.140 



18.43±0.080 
20.92±0.120 



of the best-fit model are enumerated in Table [Sj and the 
corresponding critical curves are shown in Figure [6] 

In order to scale the measured brightness back to the 
source plane and estimate the luminosity of the source, 
we need to know the average magnification of the source 
due to lensing. We estimate the magnification uncer- 
tainty through a Monte-Carlo simulation, in which we 
compute many lens models, in each one drawing a set 
of model parameters from steps in the MCMC that are 
within [Xmi„,Xmin + 2] of the best-fit model. The average 
magnification of the counter-image and its uncertainty as 
predicted by these models is 2.04 ± 0.16. The magnifi- 
cation of the giant arc is highly sensitive to its location 
relative to the critical curve in the image plane, and con- 
sequently very sensitive to the details of the model. Areas 
of the arc that lie closer to the critical curves are more 
highly magnified. The resolution of the current imaging 
data is insufficient to obtain separate magnifications for 
the 3 merged images that make up the giant arc or for 
any of the visible substructure. We therefore estimate an 
average magnification for the whole arc from the more ro- 
bust magnification estimate of the counter-image and the 
relative sizes of the arc and counter-image on the sky. We 
count the number of pixels above a brightness treshold 
in the /-band data for both; for a range of surface bright- 
ness tresholds fj, between 24.32 and 25.07 mag arcsec^-^, 
we find a relative size ratio of 8.45 ± 0.15. After mul- 
tiplying this by the magnification of the counter-image, 
we estimate the average magnification of the giant arc to 
be 17.2 ±1.4. 

5. STELLAR ENERGY DISTRIBUTIONS 

5.1. Modeling procedure 

The spectral energy distribution (SED) of a galaxy is 
governed by its star formation and metal enrichment his- 
tory in conjunction with its current stage of stellar evolu- 
tion and dust attenuation and as such contains informa- 
tion on the galaxy's stellar mass, age, star formation his- 
tory, metallicity and dust content. It is unfortunately not 
an easy task to extract these parameters with a high de- 
gree of certainty from an observed SED, a problem which 
gets increasingly worse towards higher redshifts. Recent 
studies have shown that the choice of stellar evolution 
prescription is the largest systematic uncertainty in the 
interpretation of a galaxy's SED, overshadowing the ran- 
dom errors which arise from the photometric uncertain- 
ties. It is not a goal of this paper to conduct a complete 
comparison (see e.g. Maraston 2005; Conroy & Gunn 
2009; Muzzin et al. 2009; Wuyts et al. 2009), but it is 
important to understand the main aspects. One of the 
major challenges remains the treatment of thermally pul- 
sating asymptotic giant branch (TP-AGB) stars, which 
contribute significantly to the restframe near-IR emission 
for galaxies in the age range 0.2 < t < 2 Gyr (jMaraston, 




Fig. 6. — Strong lensing interpretation shown on a composite 
color image of uBg, rlz, and JHKs- The critical curves of the 
best-fit lens model are overplotted in red: the outer curve is the 
tangential critical curve, and the inner one is the radial critical 
curve. The cuspy and oval yellow curves are the tangential and 
radial source-plane caustics, respectively. The three images that 
compose the giant arc RCSGA 032727-132609 are labeled with Aa, 
Ab and Ac. We identify the faint blue feature labeled with (A) as 
a part of the same source, but it is not included in the analysis 
presented in this paper. The counter-image is labeled with C. The 
inset shows a reconstruction of the source from the counter-image. 
Its location in the source plane is indicated by a cyan rectangle. 
We identify other likely lensed sources, and mark these candidates 
with Sl-7. We note that SI and S2 are r-dropout galaxies, S3 and 
S4 are likely a three-image system of the same source, and S7 is 
likely a radial arc. 



2005 1 . Since the near-IR traces a galaxy's old stellar pop- 
ulation, this is a key issue in the study of stellar masses. 
The updated version of the stellar evolution models by 
Chariot & Bruzual (hereafter CB07, kindly made avail- 
able by the authors) includes a revised treatment of 
TP-AGB stars and brings them into closer agreement 
with t he other widely used class of models by |Maraston| 
( 2005 1 . We use the models based on a Chabrier initial 
mass function ( Chabrier 200^ It is well established that 
a Salpeter IMF ( |Salpeter||1955| ) overpredicts the number 
of stars less massive than 1 M© , increasing the inferred 
stellar masses and star formation rates by a factor of r^ 
1.5-2 (see e.g. Papovich et al. 2001). We investigate the 
infiuence of a Salpeter IMF for RCSGA 032727-132609 
in § |6.3| to allow comparison with other studies. The 
dust extinction is governed by the Calzetti extinction 



TABLE 3 

Best-fit lens model parameters 



Mass 



RA 

(h m s ±") 



Dec 



(dog) 



(kpc) 



rcut 
(kpc) 



(km s^-'^) 



Cluster (PIEMD) 03:27:27.310 ±0.2 

L* galaxy (PIEMD) 



-13:26:23.65 



+ 0.5 r, fjr, +0.07 

-0.3 "-"^ -0.06 



±0.4 5.3 ±1.2 

[0.15] 



281 
5.2 



+ 13 
-115 
+0.7 
-3.2 



1094 

87 



+7 
-21 

+32 

-29 



Note. — The ellipticity is expressed as e — (a — b )/(a + 6 ). ^ is measured north of West. Error bars correspond to 1(T 
confidence levels as inferred from the MCMC optimization. Values in square brackets are not optimized. The location and the 
ellipticity of the matter clumps associated with the cluster galaxies were kept fixed according to the light distribution. 



3. 

Mi 



law (Calzetti et al. 2000L which has been derived for lo- 



cal starbursts. 'i'he validity of this extinction law at high 
redshift has yet to be tested adequately. It has been sug- 
gested that the Calzetti law overpredicts dust extinction 
for young LBGs fages < 100 Myr) when compar ed to 
other diagnostics JReddy et al ^[20060 I Siaiiae taljIgOOS 
2009l|Reddy et al.plO[ ). IMuzzmeTal | ( |200D[ ) iincf altf 
2U7o systematic etiect on stellar population parameters 
from the choice of dust extinction law, small compared 
to the systematic errors connected to the choice of stellar 
evolution model. 

Once the model set has been chosen, we have to worry 
about the internal degeneracies of a spectral energy dis- 
tribution. Age, metallicity and dust all tend to affect 
a galaxy's spectrum in similar ways; the SED will over- 
all be redder for a more metal rich stellar population, 
an older stellar population or more dust extinction. To 
decrease the uncertainties of the important stellar pop- 
ulation parameters such as stellar mass and age which 
originate from this degeneracy, the metallicity is often 
treated as a known parameter. We use solar metallic- 
ity as a default choice, which is also most suitable for 
comparison with the literature. Lower metallicities are 
investigated in § 5.3[ The remaining degeneracy between 
the age of a stellar population and its reddening can be 
broken if near-IR photometry is included in the SED to 
cover the age sensitive 4000A break. 



We use an update d version of the code Hyper 



pd^ 

zonella et al. 2000), kindly made available to us by 
M. Bolzonella, which performs SED fitting at a fixed 
spectroscopic redshift. The code matches the observed 
broadband SED to a set of template spectra through 
a maximum-likelihood procedure, quantifying the good- 
ness of fit via x^-statistics. It is important to note that 
the code does not interpolate on the template grids; the 
input template set must be densely populated. We use 
the solar metallicity templates from CB07 for a range 
of exponentially declining star formation histories of the 
form SFR{t) ^ exp(— i/r), with e-folding times r — 0.01, 
0.05, 0.1, 0.2, 0.5, 1 and 2Gyr, as well as continuous 
star formation models (CSF). For the main analysis, we 
do not consider more complex star formation histories 
such as multiple component models with several short- 
duration bursts, since the data do not place strong con- 
straints on even the single-component models. We only 
briefly investigate an extreme two-component model in 
the next section to determine an upper limit on the stel- 
lar mass from a maximally old underlying burst. Ages 
are allowed to vary between 0.1 Myr and the age of the 
Universe at z = 1.701 and extinction is constrained to 



heE{B-V) <0.75. 

The procedure outputs the best-fit spectral energy dis- 
tribution and the corresponding stellar population pa- 
rameters. Independent estimates are obtained for the 
arc and the counter-image and can be used as an inter- 
nal consistency check of the photometry as well as the 
SED modeling procedure. The age i, reddening factor 
E{B — V) and star formation history r depend solely on 
the galaxy colors, and should be similar for both images 
of the same background galaxy. The stellar mass and cur- 
rent star formation rate have an additional dependence 
on the lensing magnification and are corrected for the 
magnification factors obtained in § |4j The uncertainties 
related to the lens model need to be taken into account 
when comparisons are made for these two parameters. 

5.2. Results 

Figure [7| shows the best-fit SEDs for the arc and 
counter-image found from the set of models described 
above. The fits are good, with a reduced x^ of 0.87 
and 0.89 for the arc and counter-image respectively. It 
is important to supplement the best-fit stellar popula- 
tion parameters with parameter confidence intervals al- 
lowed by the photometric uncertainties. We create 1000 
fake realizations of the observed SED by perturbing each 
broadband magnitude measurement in a manner consis- 
tent with its errorbars. This set of fake SEDs is fit in 
exactly the same manner as described above for the ob- 
served SED; bad fits with x^ > 3.0 are excluded and 
the median and standard deviation are deduced for each 
stellar population parameter. The results are listed in 
Table ID The stellar mass and current star formation 
rate have been corrected for the magnification factors 
of 17.2±1.4 for the arc and 2.04±0.16 for the counter- 
image. The correspondence between the stellar popula- 
tion parameters independently inferred from the arc and 
counter-image is encouraging. 

The stellar mass is one of the most robust parameters 
that can be determined from SED modeling, least af- 
fected by intrinsic uncertainties related to the choice of 
stellar evolution model. However, there are other effects 
to be taken into account. First of all, the stellar mass 
reported from SED modeling corresponds to the mass 
turned into stars by the age of the galaxy based on its star 
formation history. This is an overprediction since during 
the galaxy's lifetime stars die and leave remnants whose 
mass is smaller than the initial one (Renzini & Ciotti] 
1993| For the relatively young age of RCSGA 032727^ 
132B09, this mass returned to the interstellar medium 
(ISM) will be a very small percentage. 

The biggest mass bias is caused by the use of single- 
component star formation history models. These at- 




Observed wavelength [jum] 



Fig. 7. — Best fit stellar energy distributions for the arc and 
counter-image. Photometry in u, B, g, r, I, z, J, H and ii'a-bands is 
overplotted as datapoints with la errorbars. 

tribute the emission at all wavelengths from rest-frame 
UV through near-IR to a single, smoothly declining or 
continuous episode of star formation, which is not a very 
realistic picture of the episodic star formation histories 
that arise in a hierarchical model of galaxy formation. 
When a young episode of star formation is superimposed 
on an older burst that peaked sometime in the past, the 
stellar mass contribution of this older stellar population 
will likely not be captured by single-component modeling 
of the overall SED. In order to set a rough upper limit 
on the range of stellar masses a llowed by the S ED of RC- 
SGA 032727- 132609, we follow [Shapleyet al.| ( |2005| ) and 
Daddi et al. ([2004b) in combinmg a very young model, 
which dominates the emission at restframe UV wave- 
lengths, with a maximally old underlying burst, which 
dominates the near-IR emission and whose age is only 
limited by the age of the Universe at z = 1.701. Specif- 
ically, we scale a young (i = 10 Myr) CSF model to the 
observed r-band magnitude, subtract this model from the 
observed SED, and scale the maximally old {t — 3 Gyr, 
T = 100 Myr) model to match the residual near-IR mag- 
nitudes. The stellar mass is the sum of the two compo- 
nents, largely dominated by the older stellar population, 
and is found to be about twice as large as our best-fit 
single-component mass (after correcting for the mass re- 
turned to the ISM by supernovae). Similarly, the age 
inferred from SED modeling represents the duration of 
the current episode of star formation and could be larger 
for an underlying older stellar population. 



5.3. Constraints on metallicity 

We have assumed solar metallicity for the SED mod- 
els to remove one parameter in the degeneracy between 
age, dust extinction and metallicity that arises from the 
spectral slope. It is of interest to explore the effect of sub- 
solar metallicity models on the best-fit stellar population 
parameters. The presence of the CIII] nebular emission 
line could point to a metallicity similar to or less than 
0.4 Zq found for cB58 ( [Pettini et al]|2000| which does 
not show this line in its spectrum (see § 2?2). We fit the 



same range of models (Chabrier IMF, Calzetti dust ex- 
tinction, T = 0.01, 0.05, 0.1, 0.2, 0.5, 1.0, 2.0 Gyr or CSF) 
with 0.4 and 0.2 Zq metallicities. We find acceptable fits 
for each metallicity, the best-fit SEDs are shown for the 
counter-image in Figure [81 

Based on x^-statistics, the 0.2 Zq models produce a 
worse fit, but the 0.4 Z© best-fit model reports a bet- 
ter x^ of 0.66 and 0.63 for the arc and counter-image 
respectively, compared to 0.87 and 0.89 for the default 
solar metallicity. The stellar population parameters from 



the 0.4 Zq models are listed in Table 



m They 



are overall 



consistent with the confidence intervals deduced earlier 
from the solar metallicity models; the smaller ranges al- 
lowed for the age and stellar mass confirm that the 0.4 Zq 
models produce overall better fits to the observed SED. 
Erb et al. (2006) have investigated the mass-metallicity 
relation for galaxies at z ~ 2 and find a metallicity of 
^ 0.5 Zq for a stellar mass of 10^" Mq, consistent with 
our results. From Figure [^ it is obvious that the SEDs 
start to diverge considerably in the mid-IR region, such 
that inclusion of IRAC data to the SED could distinguish 
between these models. This will be explored in a future 
paper. External metallicity information from spectral 
emission lines would also be very helpful. 








Observed wavelength [jum] 



Fig. 8.— Best-fit SEDs from Zq, 0.4 Z© and 0.2 Z© metallicity 
models for the counter-image. Photometry in u, B, g, r, Iz, J, H and 
A's-bands is overplotted as black datapoints with la errorbars. 



6. CONTEXT 

6.1. Luminosity 

An often quoted measure of the brightness of a galaxy 
is its comparison to L* , the luminosity of the c haracteris- 
tic bre ak in the Schechter luminosity function ( Schechter 
There have been a number of recent studies of 



1.7. Arnoutset al. (20051 



1976[ ). 
the luminosity function at z 
provide an estimate of the UV luminosity function at 
1500A for a sample of galaxies at 1.75 < z < 2.25 with a 
characteristic absolute magnitude M* = — 20.33 ± 0.50. 
This i s consistent with a recent study by |Oesch et al.| 
( 2010| which reports M* = -20.34 ± 0.29 for galaxies at 
5 < z < 2.0. We compare this to the absolute mag- 
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TABLE 4 

Stellar population parameters 



ARC 

Z© 0.4 Zq 



COUNTER-IMAGE 
Zq 0.4 Zq 



Age [Myr] 91±65 91±27 143±90 72±35 

E{B -V) < 0.035 < 0.110 < 0.025 < 0.045 

T [Myr] 10-50 10-50 10-50 10-50 

log(M/MQ) 10.07±0.08 10.05±0.07 lO.OlitO.lO 9.95±0.07 

SFi? [M0yr-i] < 35 < 77 < 38 < 41 



nitude at rest-frame 1500A of the counter-image of RC- 
SGA 032727-132609, derived from the best fit SED via 
Mi5oo = "11500 - 5 log(i:>t/10pc) -t- 2.5 log(l + z).We use 
the counter- image as the best representative of the source 
because of its more robust magnification estimate. Tak- 
ing into account the lensing magnification of 2.04 ±0.16, 
we find M1500 = —22.37 ± 0.09. Relative t o the charac- 
teristic M* = -20.34 ± 0.29 postulated by |Oesch et al. 



( 2010 ), we find RCSGA 032727-132609 at a luminosity 
oT^il.SL*. 

6.2. Stellar mass 

Similar to the comparison to L* made above, we 
can compare the stellar mass of RCSGA 032727-132609 



to M* 



the mass of the characteristic break in the 



Marchesini et 



Sch echter mass function ( Schechter 1976 ) 
al. ( 2009 1 present a thorough analysis of the evolution of 
the stellar mass function from z ~ 4 to the present. They 
devote much attention to the systematic uncertainties in- 
duced by the assumptions made in the SED modeling of 
the stellar masses. For the best-fit Schechter function de- 
rived with SED assumptions consistent with ours (CB07, 
Kroupa IMF-consistent with Chabrier, Calzetti extinc- 
tion, Zq) for the redshift range 1.3 < z < 2.0, they find 
log(M*t^,/M0) = 10.80 ± 0.05. When comparing this to 
the stellar mass of log(M/MQ) = 10.01 ± 0.10 obtained 
from the counter-image, we find RCSGA 032727-132609 
at a stellar mass of 0.16 ± 0.04 M*^^,. 

6.3. Comparison to representative galaxy samples 

In addition to the statistical measures discussed above, 
we want to place RCSGA 032727-132609 in reference to 
the overall galaxy population at this redshift that ex- 
ists in current studies, to determine to what extent the 
stellar population parameters found in this highly mag- 
nified galaxy can be generalized to the known popula- 
tion as a whole. The most representative sample for 
comparison are the UV selected star-forming galaxies 
at 1.4 < z < 2.5 collected by Steidel and collaborators 



i^delberg er et al.|2004[|Steidel et al.|2004||Shapley et al 
1^05). They have adapted the color criteria from the suc- 
cesful LEG continuum break technique at z '^ 3 (Steidel 
1996; 2003) to lower redshifts, defining 'BM' and 'BX' 
color crite ria for 1.4 < z < 2.0 and 2.0 < z < 2.5 re- 
spectively ( [Adelberger et al.|2004[ ). |Reddy et al]j2006a I 
present SED modeling tor a large sample of 212 spec- 
troscopically confirmed 'BM/BX' galaxies. Because of 
the overall faint magnitudes, the star formation histories 
remain basically unconstrained by SED modeling and 
they focus on stellar masses as a more robust diagnos- 
tic. Specifically for the subsample of 51 'BM' galaxies 
at z = 1.72 ± 0.34, they find an average Salpeter stellar 



mass of (log(M/Mo)) = 9.86 ±0.57. Their sample spans 
a large dynamic range in age (few Myr to 2.5 Gyr), red- 
dening {E{B ~ V) — 0.0 — 0.5) and star formation rate 
(5fi?=0-914Moyr-i). 

When fitting the SED of the counter-image of RCSGA 
032727-132609 with solar metalhcity Salpeter IMF mod- 
els, we find a stellar mass of log(M/M0) = 10.26±0.10, a 
factor of 1.8 higher than what was stated above based on 
the Chabrier IMF, in agreement with the known overpre- 
diction of the number of low mass stars by the Salpeter 
IMF. This is consistent with the average stellar mass 
found by |Reddy et al. (2006a I. The best-fit age, redden- 
ing and star formation rate of RCSGA 032727-132609 
also fall into the large ranges found for this comparison 
sample. 

7. SUMMARY 

This paper reports on the discovery of an exceptionally 
bright and extended star- forming galaxy at z = 1.701, 
strongly-lensed by a foreground cluster discovered in the 
RCS2 survey. The giant arc is ~ 3.5 times brighter than 
cB58 and extends over 38 "on the sky. We measure a ve- 
locity dispersion of 988±122kni s~^ for the cluster and 
estimate a virial mass of M200 ~ 1.1 x IO^^MqIi^q, 
in accord with the large Einstein-radius of ~ 17'.'8. 
A lens model is constructed for the cluster using the 



publicly-available software LENSTOOL (Julio et al. 
magnification factor of 2.04±0 






20071 
6 tor 



and results in a 
the counter-image and an estimate of 17.2±1.4 for the 
average magnification of the giant arc, based on the rel- 
ative sizes of the arc and the counter-image on the sky. 
Higher resolution imaging is required to create a robust 
lens model for the giant arc which correctly explains the 
apparent substructure. 

Careful measurements of consistent magnitudes for 9 
bands of photometry ranging from u to Kg produce a 
well-constrained spectral energy distribution. SED fit- 
ting is based on CB07 models with a Chabrier IMF 
and Calzetti dust extinction. The systematic uncer- 
tainties which arise from the choice of stellar evolution 
model, IMF and extinction law can be large and of- 
ten exceed the statistical uncertainties from photomet- 
ric errors. Our best-fit model places RCSGA 032727- 
132609 at a metalhcity of 0.4 Z© with a moderately young 
age, t — 80 zt 40 Myr and a relatively small amount 
of dust extinction, E{B ~ V) < 0.11. Taking into ac- 
count the lensing magnification, we find a stellar mass 
of M* ~ 10^° M0 and a current star formation rate 
SFR < 77M0yr~^. The agreement between the stellar 
population parameters independently inferred from the 
giant arc and counter-image, both images of the same 
background galaxy, is a good consistency check of the 
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photometry, SED fitting procedure and lens modeling. 
Single-component SED models report ages and stellar 
mass measurements for the current episode of star for- 
mation. When allowing for an episodic star formation 
history, twice this mass could be present in an older un- 
derlying burst without affecting the SED detectably. 

A comparison to the characteristic L* and Ml^.^J. of 
the Schechter luminosity and mass functions at z ^ 
1.7 places RCSGA 032727-132609 at a luminosity of 
6.5 ± 1.8 L* and a stellar mass of 0.16 ± 0.04 M*^^,, sug- 
gesting a low-mass starbursting galaxy. The Salpeter 
stehar mass estimate of RCSGA 032727-132609 is con- 
sistent with the average stell ar mass for a sample of 'BM' 
galaxies at z = 1.72 ± 0.34 ( [Reddy et ar]|2006a[ ) and its 
age, reddening and current star formation rate generally 



fall into the broad range of stellar population parameters 
found for this comparison sample. The large magnifica- 
tion and spatial extent of RCSGA 032727-132609 provide 
a unique opportunity to spatially resolve one representa- 
tive example of the variety of stellar populations found 
at these redshifts and study it in great detail. This paper 
has presented only our first level of understanding of this 
unique galaxy. 
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